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Abstract 

In this paper we solve one dimensional SU(3) bosons with repulsive (^-function 
interaction by means of Bethe ansatz method. The features of ground state 
and low-lying excited states are studied by both numerical and analytic meth- 
ods. We show that the ground state is a SU(3) color ferromagnetic state. The 
configurations of quantum numbers for the ground state are given explicitly. 
For finite N system the spectra of low-lying excitations and the dispersion re- 
lations of four possible elementary particles (holon, antiholon, cr-coloron and 
o;-coloron) are obtained by solving Bethe-ansatz equation numerically. The 
thermodynamic equilibrium of the system at finite temperature is studied 
by using the strategy of thermodynamic Bethe ansatz, a revised Gaudin- 
Takahashi equation which is useful for numerical method are given . The 
thermodynamic quantities, such as specific heat, are obtain for some special 
cases. We find that the magnetic property of the model in high temperature 
regime is dominated by Curie's law: x«l /T and the system has Fermi-liquid 
like specific heat in the strong coupling limit at low temperature. 
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I. INTRODUCTION 



One of the main goals of theoretical physics during the past 40 years is to understand 
quantum systems involving many particles. Once the interaction between those particles 
were taken into account, the problem becomes complicated. Meanwhile, as long as their in- 
teraction is not sufficient weak, the perturbative methods that were powerful in many quan- 
tum mechanics text book become unreliable. In one dimension, various non-perturbative 
methods were proposed, among which the impact of exactly solvable theoretical models are 
undeniable. Particles with 5-function interaction is a simple but interesting model. Lieb and 
Liniger [1] first solved a Bose system under periodic boundary condition in the case of spin-0 
or in the absence of internal degree of freedom. The method they used is nowadays referred 
as coordinate Bethe ansatz. The extension of periodic boundary condition to the boundary 
condition of potential well of infinite depth were made by Gaudin [2] and Woynarivich [3] . 
The first attempt to develop the method applied in [1] to deal with spin-1/2 fermions was 
made by McGuire [4] who can, however, deal with the case of only one spin-down and the 
other spins keep up. Further step of considering two spin-down and others spin-up was made 
by Flicker and Lieb [5]. Gaudin [6] and Yang [7] successfully solved the problem for arbitrary 
number of spins in the state of spin-down. In Yang's paper [7], the first non-trivial case of 
Yang-Baxter equation was introduced. Actually, the strategy for general multi-component 
systems was proposed in [7] though the explicit solution was give only for spin-1/2 particles 
(it is Fermi system then). 

The literature in [7] was later extended by Sutherland [8] to the called any irreducible 
representation of permutation group. Actually, both Yang and Sutherland adopted antisym- 
metric wave function under permutating indistinguishable particles. Thus Yang solved the 
problem of 2-component fermions and Sutherland solved N-component fermions by means 
of coordinate Bethe ansatz. As 2-component system is mostly associated with "spin-1/2" 
system which is conventionally referred to Fermi system, the coordinate Bethe ansatz has 
not been employed to 2-component Bose systems till recently [10]. 
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Along with the developments of quantum inverse scattering method, Kulish [11] discussed 
multi-component nonlinear Schrodinger equation in terms of quantum inverse scattering 
method (QISM) in order to re-derive the Bethe-ansatz equations of Yang [7] and Sutherland 
[8]. The author explicitly formulated 2-component case and conjectured that Sutherland's 
results would be obtained by repeating his procedure n — 2 times. Actually, it is not possible 
because the author [11] adopted commutation (instead of anti-commutation) relations but 
the system both Yang and Sutherland considered is Fermi system. It is now clear that the 
first quantization form of the system which Kulish considered ought to be a system of SU(n) 
bosons with 5-function interaction. The QISM was also employed to nonlinear Schrodinger 
equation for graded matrix but it breaks the Yang-Baxter relation at first [12] that was 
noticed and overcome later on [13]. 

Although the Bethe-ansatz equations for 2-component bosons was early formulated by 
Kulish, the nature of ground state and the property of low-lying excitations has never been 
exposed till the paper of Li et. al. [10]. We know that not only the 2-component Bose 
gas can be formed in magnetically trapped 87 Rb [14], but also a 3-component Bose gas can 
be produced in an optically trapped 23 Na [15], it will be valuable to study the model of 
3-component Bose system. In present paper, we study a system of three-component bosons 
with SU(3) symmetry in one dimension. On the basis of Bethe-ansatz equations, we discuss 
the ground state, low-lying excited states and the thermodynamics of the system at finite 
temperature, and try to obtain thermal coefficients for some special cases. Our paper is 
organized as follows: In the following section we introduce the model and the corresponding 
Bethe-ansatz equations for charge rapidity and color rapidities. In Sec. Ill, we explicitly 
show that the ground state is a color ferromagnetic state and how the quantum numbers 
in Bethe-ansatz equations should be taken for the ground state. In Sec. IV, We study 
the low-lying excited states extensively by analyzing the possible variations in the sequence 
of quantum numbers. Numerical results of energy momentum spectra for each excitation 
are given. Furthermore, the dispersion relation of four possible elementary particles are 
obtained. In Sec. V we discuss the general thermodynamics of the system with the strategy 
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of thermodynamic Bethe ansatz (TBA) which was proposed by C.N. Yang and CP. Yang 
when they study Bose gas with 5-function interaction in one-dimension [20]. In Sec. VI we 
discuss the system for some special cases and obtain some analytical results. In Sec. VII a 
brief summary is given. 



II. THE MODEL AND ITS BETHE-ANSATZ SOLUTION 

We consider interacting SU(3) Bose field in a one dimensional ring of length L. The 
model Hamiltonian of the system reads, 



7io = J dx 



a A a,b 



(1) 



where natural unit is adopted for simplicity. Here c is the coupling constant and a,b = 
1,2,3 (we call colors hereafter) denotes the three states that carry out the fundamental 
representation of SU(3) group. The fields obey the following commutation relations, 

{r a (x),My)} = E S abS(x -y-nL). (2) 

n 

In the terminology of group theory, the three states |1), |2) and |3) are labeled by weight 
vectors, (1/2,0), (—1/2,1/2) and (0,-1/2) respectively. The two su(2) subalgebra in the 
s«(3) Lie algebra are [T+T"] = 2T Z and [U + , U~] = 2U Z . With the help of those "flipping" 
operators, U ± and T ± , we can generate the three states from the highest weight state |1), 
i.e., 

T"|l) = |2>, U~\2) = |3>. 
With additional commutation relations defined by 

V+=[T + ,U + ], V- = [T-,U-], 

the Chevalley bases of the su(3) Lie algebra consists of eight generators 
{T ± ,U ± ,V ± ,T Z ,U Z }. 
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In the domain with Xi ^ Xj, the Hamiltonian (1) reduces to the one for free bosons 
and its eigenfunction is therefore just the superposition of plane waves. When two particles 
collide with each other, a scattering process occurs. The coordinate Bethe ansatz embodies 
that this process is purely elastic, i.e., exchange of their momenta. So for a given momentum 
k = (ki, k<ii ■ ■ ■ , /cat), the scattering momenta include all permutations of the components of 
k. Thus for the case of N bosons, due to the Hamiltonian is invariant under the action of 
the permutation group Sn, one can adopt the following Bethe-ansatz wave function, 

* a (x) = £ A a (P,Q)e^ x \ x e C(Q), (3) 
Pes N 

where a = (ai, a 2 , . . . , a at), Q>j denotes the color label of the jth particle; Pk stands for the 
image of a given k := (ki, k 2 , ■ ■ • , k N ) by a mapping P G Sn and the coefficients A(P, Q) are 
functional of P and Q where the Q denotes a permutation of the coordinates which defines 
a region with < xq 1 < xq 2 < . . . < xq n < L. For a Bose system, the wave function is 
supposed to be symmetric under any permutation of both coordinates and color indices, i.e. 

(WV) a (x) = * a (x), (4) 

where IP : {a±, cij, ...} i— > {ai, %, ...} and (rP\I/) a is well defined by 
^/ nJa (Wx). Furthermore, using the identity {Pk\Yi % Qx) = (U l Pk\Qx) and rearrangement 
theorem of group theory, we have the following consequence from (4): 

A a (P,WQ) = A U]a (WP,Q). (5) 

The ^-function term in the Hamiltonian (1) contributes a boundary condition across the 
hyper-plane x Qj = x Qj+1 , 

i{{Pk) 3 - (Pk) j+1 )[A a (P, WQ) - A a (WP, WQ) - A a (P, Q) + A a (WP, Q)} = 2c[A a (P, Q) + A a (WP, Q)}, 

(6) 

By making use of the relations (5) and (6) together with the continuity condition, we can 
obtain the following relation 
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where P J permutates the color labels of bosons located at x Qj and x Qj+1 . 

Applying the periodic boundary condition • • , xq v • • •) — • • , Xq j + L, ■ ■ •) and 
making use of the standard procedure of quantum inverse scattering method [8,11,16], one 
can obtain the Bethe-ansatz equations 

jk iL _ "kj-k + ic™ kj — \ v — ic/2 



e 



— TT — — — — — TT 

1=1 — ki — ic ^ kj - \ u + ic/2 ' 
-iX- A 7 — ki — ic/2 it A 7 — + ic 
2=1 A 7 - ^ + ic/2 ^ A 7 - \ v - ic 

W " A- /,,, ic/2 



X 



n 



t=i A 7 -/i a + ic/2' 
1 = _ yr up - K-ic/2 ¥1 Hp- H a + ic 
i=i Hp- K + ic/2 ^ fj,p - fi a - ic 

The A and ji are SU(3) color rapidities. There are M — M' particles in the state |2), 
M' in |3) and N — M in |1). We would like to mention here that the state obtained 
above is the highest weight state among the multiplet of SU(3) representation labeled by 
(N/2 + M'/2 — M, M/2 — M'). The other states in the multiplets can be generated by iterate 
application of the flipping operators T~ and U~ . 

Taking logarithm of Eqs. (8) we have secular equations, 

N M 

kjL = 2nlj + i( A; i - k i) + E @-i/2 (^i - K), 

1=1 v=l 
N M 

2ir J 7 = £ 0_i /2 (A 7 -h) + J2 ©i(A 7 - A,) 

1=1 v=l 
M' 

+ J2 @-l/2 (A 7 
a=l 

M M' 

2nJ'p = ]T 0-i/2(^ -K) + Y1 0i(^ - /U (9) 

v=l a=l 

where Q n (x) = — 2 tan _1 (:r/nc). The quantum number Ij for charge rapidity kj takes integer 
or half-integer depending on whether N — M is odd or even. The quantum number J 7 and 
J'p for SU(3) color rapidities A 7 and /ip take integer or half-integer depending on whether 
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iV — M — M' and M — M' is odd or even respectively. Once all roots {kj, A 7 , pp} are solved 
from the above equations (9) for a given set of quantum numbers {Ij, J 7 , J'p}, the energy 
and momentum will be calculated by 

N 

' N M M' 

EVEVEA , (io) 

J'=l 7=1 13=1 



2tt 



where the second equation of Eqs. (10) is obtained from Eqs. (9) directly. 

For a state with real roots (k, A,//), we may define the distribution densities p(k), er(A) 
and uj{p) 



p{kj) 




kj), 


cr(A 7 ) 




- A 7 ), 






- A*/?)- 



(11) 

In terms of those densities the energy and momentum become 

£/L = /" k 2 p{k)dk, 

p/L = J kp(k)dk. (12) 

While JV, M and M' are determined by 

N/L = J p(k)dk, 
M/L = J a(X)d\, 

M'/L = J uj(p)dp. (13) 

As the SU(3) "magnetic" field is characterized by two parameters Hi and H 2 , the Zeemann 
term is given by 

H zee = Hx(N -2M + M')/2 + H 2 (M - 2M')/2 



1 p{k)dk + V 2 J / a(X)dX 



ui{p)dp. 



[H, - 2H 2 )L 



(14) 
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III. THE GROUND STATE 



It is easy to show that the first equation of Eqs. (9) is a monotonously increasing function 
of kj, that is if ij < Ij we have fcj < kj. So the configuration of {Ij} for the ground state is 
given by succussive integers or half integers that is symmetrically arranged around zero, i.e., 
Ij + i — Ij = 1. In order to observe the properties of { J 7 , J'p}, it is useful to investigate the 
Eqs. (9) in the weak coupling limit c — > 0. Due to 0± n (x) — * =)=7rsgn(a;), Eqs. (9) become 

V M 

21 j = kjL/7r + sgn(/cj ~ki) -J2 s g n (^j - 

N M 

1=1 v=\ 
M' 

+ J2 sgn(A 7 -// a ), 

a=l 

M M' 

2 J'p = sgn(/i /3 - A„) - s S n (/^ - /V>- ( 15 ) 

i/=l a=l 

We can choose the subscripts of the rapidities kj, A 7 , /ip in such a way that Ij, J 7 , J'^ are all 
ranged in an increasing order. Then we have 

L M 
2(I j+1 -Ij - 1) = -{kj+i - kj) - ^[sgn(%+i - A„) - sgn(A; j - A„)], 

71 u=l 

N M' 

2(J 7+1 -J 7 + 1) = ^][sgn(A 7+ i ~ k l) ~ S g n (^7 ~ k l)] + H[ S g n ( A 7+l ~ Z 1 ") ^ S g n (-^7 -J"a)], 
Z=l a=l 

M 

2(J'/3+i- A + !) = X)[ s g n W+i - A„) - sgn^ - A„)]. (16) 

u=l 

Therefore, if J'p+i — J'p = m, there must exist m + 1 solutions of A„ satisfying //g < 
A^ < fJ-p+i", and if J 7+ i — J 7 = n, there must be n + 1 solutions of k\ and /i a satisfying 
A 7 < ki,/i a < A 7+1 . So the existence of a Aj, between two /i's has a positive contribution 
to the density of \x (11), and vice versa for /j, to A. However, from the first equation of 
Eqs. (16), for Ij + i — Ij = n, there will be kj + i — kj = 2nn / L if there is a A 7 such that 
kj < A 7 < kj + i, otherwise kj + i — kj = 2(n — l)n/L. So a rapidity of A 7 always repels the k 
rapidity away from that value. As a result, an existing A 7 will suppress the density of state 
in /c-space at the point k — A 7 . The weaker the coupling, the more magnificent the effect 



will be (see Fig. 2). And for a given set {/}, the more A rapidities there are, the higher the 
energy is. 

It is also useful to observe Eq. (16) in the strong coupling limit. We consider two cases: 
M = and M — 1. For M — 0, the secular equation becomes 

N 

kjL = litlj + ]T Giikj - h). (17) 
i=i 

and for M — 1 we have 

N 

k)L = 2tt/; + £ - k[) + e_ 1/2 (^ - a x ). (is) 
i=i 

Here Jj — = 1/2 due to M changing from zero to one. As c — > oo, we have tan _1 (a;/c) ~ 
x/c. So the above two equations become 

2 AM 



(kj+i — kj)L 
(k' j+1 - k'j)L 



1+ Lc 
1 + 2(iV-l) 



2tt, 

2tt. (19) 



Lc 

Whence the distribution is almost a histogram. Referring to Eq. (11) the value of the density 
distribution for M = is larger than that for M — 1, which makes the Fermi momentum 
for later case to be larger than the former case so that to keep the total number of particles 
being the same. Therefore the state of M = has lower energy. 

Differing from the SU(3) fermionic model [8] and a toy model of quark cluster [9], the 
ground state of SU(3) bosonic model is no more a color singlet but a color ferromagnetic 
state. The difference is due to the distinct permutation symmetries. For iV particles, the 
ground state is characterized by a one-row iV-column Young tableau [N] whose quantum- 
number configurations are 

{-(AT - l)/2, • • • , (iV - l)/2} 
M= M' = 0. (20) 

The density of states for the ground state is plotted in Fig. 1 for various coupling with 
L = N = 41. 
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In the thermodynamic limit, the density corresponding to the configuration of quantum 
numbers of the ground state satisfies the integral equation 

Po(k) = t F K 2 (k - k')p (k')dk'. (21) 

Z7T J-kp 

Here po(k), kp are respectively the density and integration limit for the ground state, and 

*.(*) = ±- " c/2 



7r n 2 c 2 /4 + x 2 ' 
The concentration is given by 



D = N/L = [ kF p (k)dk. (22) 

J — kp 



From Eq. (21) and Eq. (22), we can determine poik) and kp. Here kp is a quasi-Fermi 
momentum because the wave function vanishes for any kj = k\ (j ^ I) as long as c ^ even 
in Bose system which can be seen from Eq.(7). The energy can be calculated by 

E /L = f kF k 2 p (k)dk. (23) 

J— kp 

1 4 

which is explicitly -ir 2 D 3 (l — D) in the strong coupling limit. In the general case one needs 
3 c 

to solve the equations numerically. We show the ground state energy for particle densities 
D = 1.0,0.75,0.5 in Fig. 3. 



IV. LOW-LYING EXCITED STATES 

The low-lying excited states are obtained by varying the configuration {Ij, J 7 , J'p} from 
that of the ground state. 

Holon-antiholon excitation. The simplest case is to remove one of / from the config- 
uration of the ground state and add a new one outside the original sequence, i.e., 

{/,} = {-(N - l)/2, m - l, m + l, (N - l)/2, /„}, 

where \I n \ > (N — l)/2 and M = M' = 0. We call this holon-antiholon excitation which 
consists of a "ho Ion" created under Fermi surface and an "antiholon" created outside it. 
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In Fig. 4, we plot the numerical results of energy-momentum spectrum for a system with 
L = N = 41 (the other part is just the mirror image of the plotted part corresponding to 
the state with p — > —p coming from the negative I n ). From the figure, we notice that there 
is a minimum in the excitation energy at p = 2n due to the fact that both 7° replaced by 
I n = (JV + l)/2 and 1% share the same energy, their momenta difference, however, is 2n. The 
overall structure of the spectrum is not changed obviously between c = 1 and c = 10. For 
a system of finite size, the gap of holon-antiholon excitations opens. In the thermodynamic 
limit, however, it vanishes. 

In the configuration of quantum numbers for the ground state (20), replacing 1% = 
(N — l)/2 by 1% = (N — l)/2 + n, n = 1,2, •• • and keeping the others unchanged, we 
obtain the dispersion relation of antiholon (Fig. 7) by solving the Bethe ansatz equations 
(9) numerically. In a similar way, replacing 7°, n — 1, • • • , TV of {I®} in turn by (N + l)/2, 
we have the dispersion relation of holon which is shown in Fig. 8. 

In the thermodynamic limit, it is plausible to calculate the excitation energy by making 
p{k) = po(k) + pi{k)/L where po(^) is the density of the ground state. By creating a hole 
inside the quasi Fermi sea k G [— kp, kp] and an additional k p > kp outside it, we have the 

Pl (k) + S(k -k)= [ kF dk' Pl {k')K 2 {k - k') 

+K 2 {k-k p ). (24) 

The excitation energy consists of two terms AE = J k 2 pi(k)dk + k 2 v = Eh{k) + s a (k p ). 
The holon energy e h and antiholon energy e a {k p ) = —Sh(k p ) are given by, 

e h {k) = -P + f kF k 2 p\{k~k)dk, 

J —kp 

p^(k,k) = -K 2 (k-k) 

/k F 
K 2 {k-k')p\{k\k)dk' . (25) 
-kp 

Holon-coloron excitation: Excitations related to the color sector is characterized by 
adding A and p rapidities into the system. The simplest excitation of this type is obtained 
by considering M = 1 which is labeled by (N/2 — 1,1/2). Comparing to the ground state, 
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the quantum number changes from half-integer to integer or vice versa. We call this type 
excitation a-coloron which is regarded as an elementary quasi-particle of the present model. 
It's quantum number takes 

Il = - N /2 + 5 1Jl (l<h<N + l), 

I J • ' 1 • (j=2,...,N), 

while Ji = Ii + m (m = 1,2, ...N — 1) so that ii < J\ < I N . This produces a N — 1 
multiplets. The excitation spectrum are plotted in Fig. 5 for a system of N = L = 21 with 
c = 10, 1 respectively. 

Adding an additional A rapidity to the color ferromagnetic ground state brings about a 
hole in the fc-sector. Now we have two-parameter excitation AE = J k 2 pi{k)dk where the 
Pi(k) solves 

Pl (k) + 5(k -k)= [ kF K 2 (k - k') Pl (k')dk' 

-K^k-X). (26) 

The energy of the holon-coloron excitation consists of two terms AE = Sh{k) + e c (X). 
The Eh is determined by Eqs. (25) and the e c defined by £ C (A) = / k 2 p\{k, A) with 

p c 1 (k,X) = -K 1 {k-X) 

+ j kF K 2 (k - k')p\{k', X)dk'. (27) 

J— kp 

The Shik) and £ C (A) are energies of holon and cr-coloron whose dispersons are shown in Fig. 
8 and Fig. 9 respectively. 

Furthermore, the overall structure for the case of c = 1 in Fig. 5 differs from the case 
of c = 10. We interpret the phenomenon as being due to the fact that the dependence of 
the dispersion relation of holon and <r-coloron on the coupling constant are different. This 
feature can be concluded from Fig. 8 and Fig. 9. When c decreases, the Eh(p) decreases 
while e c (p) increases. 
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The (j-type coloron-coloron excitation: Leaving the configuration of the ground 
state {Ij} unchanged and changing M from zero to M — 2, which corresponds to (N/2—2, 1), 
a two parameters excitation in A-sector is characterized by 

—(N - l)/2 < Jj < J 2 < (N - l)/2. 

There are totally iV(iV — 3)/2 possible choices for such type of excitation. In Fig. 6, we plot 
the numerical result of energy- momentum spectrum for a system with N = L = 21. 

The excitation energy of coloron(cx type)-coloron (a type) can be calculated by AE = 
J k 2 pl(k, Ai, \2}dk, where p1(k, Ai, A2) is determined by 

p\(k, A ls A 2 ) = -K-i(k - Ax) - - A 2 ) 

+ AT 2 (fc - k')p\{k\ Ai, A 2 )dA/. (28) 

It consists of two terms = £ c (Ai) +£ C (A 2 ), where the coloron energy £ C (A) has been give 
in the passage before Eq. (27). 

Dispersion relation of w-coloron: The fourth possible excitation involves both the 
additional quantum number J and J'. For M = 2 and M' — 1 there is no range for J' 
varying, but for large M the excitation is no more low-lying excitations. So we only show 
its dispersion relation which is described by the following configuration 

{/i}={-(JV-l)/2,...,(N-l)/2}, 
{J 7 }={-M/2,...,(M-2)/2}, 
J\=-M/2 + l,---,M/2-l, (29) 

for a given M. We plotted the dispersion relation of cj-coloron in Fig. 10 by varying J\ for 
a system with L = N = 40. The figure has a minimum around p = n when M = N/2. 

Up to now, we have discussed three low-lying excitation energies and the dispersion 
relations of four possible elementary particles, holon, antiholon, a-coloron and cu-coloron. 
We found those low-lying excitations are gapless in thermodynamic limit. 
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V. THERMODYNAMICS AT FINITE TEMPERATURE 



For the ground state (i.e. at zero temperature), the charge rapidity ks are real roots of 
the Bethe-ansatz equations (9). For the excited state, however, the A and \i rapidity can 
be complex roots [17,18] which always from a "bound state" with several other As. This 
arises from the consistency of both hand sides of the Bethe-ansatz equations [19] in the limit 
L — > oo, N — > oo. The n-string rapidity is defined by 



K J = K + (n + 1 - 2j)iu + 0(exp(-«51V)), 
K j = AC + (n + 1 - 2j)iu + 0(exp(-6N)) 

where u = c/2, j — 1, 2 • • -n. The total number of A and \i are determined by 

oo oo 

M = j>M n ; W = 

n=l n=l 



(30) 



(31) 



where M n and denote the number of A n-strings and /i n-string respectively. The Eqs. 
(9) become 

kjL = 2kIj + £ ©l(^ - *l) + E - K) 



2ttJ' 



EM^."-*i)+ E ^©^(a^-aD + e^g^/^a:-^), 
J2B nlt e t/2 (^ a -x l b )+ E A nH e_ t/2 (^-^). 



(32) 



bit 



where 



1, for t = n + I, \n — l\, 

2, for * = n + I - 2, • • • , \n - l\ + 2, 
0, otherwise. 



and 



Brdt = < 



1, for £ = n + I — 1, n + I — 3, \n — l\ + 1, 
0, otherwise. 

and the quantum numbers {Ij, J™, J'™} label the state which is no more the ground state. 
Replacing kj, A", /z™ in Eqs. (32) by continuous variables k, A, ji but keeping the summation 
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still over the solutions of these roots, we can consider the quantum number Ij, J™, J' n a as 
functions I(k),J n (X) and J' n (p) given by Eqs. (32). Take /(A) as an example, when I(k) 
passes through one of the quantum number Ij, the corresponding k is equal to one of the roots 
kj, so is for J n (X) and J' n {p). However, there may exist some integers or half-integers for 
which the corresponding k(X, p) is not in the set of roots. Such a situation is conventionally 
referred as a "hole". In the thermodynamic limit, we may introduce the densities of real k, 
A n-string and p n-string 

p{k) +p h {k) = (l/L)dI(k)/dk, 
a n (X) + a*(X) = (l/L)dJ n (X)/dX, 

luM + u h M = (l/L)dJ' n (p)/dp. (33) 
Then Eqs. (32) give rise to the following coupled integral equations, 

p + p h = ^ + J K 2 (k- k')p(k')dk' 

- W K n (k - X)a n (X)dX, 

o h n = J K n (X - k)p(k)dk 

~Y A riit ( K t {X - X')ai{X')dX' 
u J 

+ Y B mt / K t(X - fj)ui(fj)dfj,, 
it J 

"n = Y, B mt I K t (p-X)a l (X)dX 
it J 

- Y Anit I Ktifj, - p'Uip'W- (34) 

It J 

The <t„ and u n arising from the definition (34) that occur in the left hand side have been 
moved to right hand side by including the t — term in the summation. In terms of densities 
defined above, the total numbers of A and p are given by 

M/L = ]Tn f a n (X)dX, 
'/L = Y n / u n (p)dp. (35) 



M 
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In the presence of SU(3) magnetic field Hi and H 2 , we can define two type of "magne- 
tization" whose ^-components are 

T z /L = \ [ p(k)dk - / * n (X)d\ 
Z J n J 

^ n J 

U Z /L = W n I °n(\)d\-Y, n I WnWfJl. (36) 
^ n J n J 

Hence the energy contributed by the Zeemann term (14) is 

E Zee = H{F Z + H 2 U Z . (37) 
For given p(k), p h (k), cr n (\), cr^(A), u n (p) and u;£(/i) the entropy has the form [20] 
S/L = J[{p + p h ) ln(p + p h ) - p In p - p h In p h ]dk 

+ E / [K + °n) M?n + ° h n) ~ °n In O n - ^ In O^dX 
n J 

+ ^ f[(u n + w*) \n(u n + u h n ) - u n In u n - uj h n In uj^dp. (38) 

n J 

where the Boltzmann constant is set to unit. 

At finite temperature, the thermal equilibrium is obtained by minimizing the free energy 
F = E — E Zee — TS — pN where p is the chemical potential and S is the entropy of the 
system. Making use of the relations derived from Eqs. (34) 

Sp h = -5p+ f K 2 {k - k')5pdk' - ^ f K n (k - \)5a n d\, 

5a h n = f K n (X - k)5pdk - ]T A nlt [ K t (X - A')^(A')dA' + £ B nlt f K t (X - p)5uj l (p)dp, 
J it J it J 

8u h n = Y, Brat f K t (p - X)5a l (X)dX - £ A nlt f K t (p - p'^p'W , (39) 
it J it J 

and define 

^ = «(*) = e«W ^ = r] n (X) = e^/r ^ = A n (p) = e^» T . (40) 
p(k) a n (X) u n (p) 

We obtain the following conditions from the minimum condition SF = 0, namely 
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e(k) = k 2 -fi- HJ2 -T f K 2 {k- k') ln(l + e-< k '^ T )dk' - T ]T f K n (k - A) ln[l + e~^ T ]d\ 
C„(A) = n(2^i - H 2 )/2 + T f K n (X - k) ln[l + e -^/ T ]dA; + T £ A ntt / X t (A - A') ln[l + e'^'^dX' 

it J 

Uv)=n(2H 2 -H 1 )/2-Tj2B nlt f K t (fi-X) \n[l + e'^jdX 

it J 

+T £ A* f K t {ii - & ln[l + e-^' T W (41) 

A more useful version of Eqs. (41) is the recursive scheme which is a revised version of 
Gaudin-Takahashi equations, as obtained by Fourier transform. 

Tin re = k 2 - fi- HJ2 - TK 2 (k) * ln[l + k' 1 ] - if n (jfe) * ln[l + t?" 1 ], 

In m = -^sech(7rA/2 M ) * ln[(l + re" 1 )^ + %)/(! + A r 1 )]> 

4-u 

lnr/ n = ^sech(7rA/2 M ) * ln[(l + t 7ti _ 1 )(1 + ^ t+ i)/(l + A" 1 )], 
4m 

In Ai = ^sech(7rA/2 M ) * ln[(l + A 2 )/(l + j^ 1 )], 

In A n = i-sech(vrA/2 M ) * ln[(l + A n _0(l + A n+1 )/(1 + rj' 1 )]. (42) 

where * denotes a convolution. And these equations are complete by the asymptotic condi- 
tions 

Urn [ln Vn /n] = (2H 1 -H 2 )/2T, 

lim [In A n /n] = (2H 2 - H 1 )/2T. (43) 

Finally, we obtain the Helmholtz free energy F = E — TS: 

LT f 

F = jiN — — — / ln[l + e" e ]rfA;, (44) 

27T 7 

and the pressure of the system 

p =-aL-rJ ,nll+e " ]dk - (45) 

which is formally the same as Yang and Yang's expression but the equation which e fulfills 
is different. 
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VI. SPECIAL CASES 



In general, the free energy can be calculated by using formula (44), where e(k) and ( n (\) 
are determined from Eqs. (41) which can be solved by iteration. In the following we will 
consider some special cases respectively as explicit results are obtainable in those cases. 

A. Zero temperature limit 

The state at zero temperature is the ground state. When T — ► 0, the first equation of 
Eqs. (41) becomes 

e(fc) = k 2 - li - HJ2 + J K 2 {k - k')e(k')dk' 

+ W K n (k - \)( n (\)d\ (46) 

n J 

Then the Fermi surface is determined by e(fcp) = 0. Since there is no hole under Fermi 
surface, we can take the ratio k = p h /p as zero when k G [—ftp, ftp]. As a result, it is 
easy to see from Eqs. (42) that r\ n = A n — > oo. That is M = M' = 0, the state is color 
ferromagnetic state. This is consistent with the conclusion obtained in Sec. IV. Then Eq. 
(46) can be rewritten as 

e (jfe) = k 2 - p - HJ2 + j kF K 2 (k - k')e (k')dk', (47) 

which gives the solution of dressed energy [21], and the ground-state energy can be given in 
terms of eo 

E /L = J- f kF e (k)dk. (48) 
whose dependence on the coupling constant is shown in Fig. (10). 

B. High temperature limit 

In the high temperature limit T — > oo, we can assume that all functions rj n (X) and A n (/i) 
are independent of their corresponding parameter. Due to lim„^o i^ sec Mf^) = E^ 8 - 
(42) can be written as follows, 
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i£ = (1 + 0(1 + *»)/(! + 

rj 2 n = (l + r jn . 1 )(l + r jn+1 )/(l + A^ 1 ), 

A? = (l + A 2 )/(l + 0, 

A* = (1 + A n _ x )(l + A n+1 )/(1 + O- (49) 

with the asymptotic conditions Eqs. (43). 

Performing the Fourier transform to Eqs. (34), we get the solution of the densities of A 
n-strings, 

0i + a\ — — sech[7rA/2-u] * [p + o\ + o^], 
4-u 

^ + ^ = ^Sech^^/ 2 ^] * + <?n + l + "n\, 

uj\ + Jl = -3-sech[7rA/2-u] * [ujL 1 + aA 

uJ n + uJ^ = ^sech[7rA/2 M ] * [uj^ + uj* +1 + a n }. (50) 

If we assume that a n , and u n ,u^ are independent of A and « respectively, or let u — 0, 
we have the following relation 

rt Z Z ,„ Z 

E ™ B = 1 £ ™« - nj ^ 1 ^e n ^ !T (51) 

n z n Z 

where n m and n m i are the maximal length of A string and « string respectively, and Qi = 
2Hi — H 2 , fi 2 = 2H 2 — Hi. In the absence of external field Hi and H 2 , it is easy to obtain 
M' = M/2 = N/3 which means there are N/3 particles in each internal state. Then we 
can also infer that the contribution of the internal degree of freedom to the entropy per site 
must be S = In 3 which follows from the fact that the internal degree of freedom per particle 
is three. 

If the external field is small, expanding Eq. (51) for small field and integrating the 
equation over A and «, we get the SU(3) magnetization of the model, 



T z M m 



L 2L 



n m tti I f n m VLi \ 2 
T 2 V T ' 
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El 

T 



M' 
2L 



n m >Q 2 1 f n m >VL2 
T 2 V T 1 



(52) 



where M m and M m / are the total number of rapidities in A n m -strings and fi n m /-strings re- 
spectively. The first term in the parentheses of both equations arises from self- magnetization, 
while the others are contributed by external field. Eq. (52) indicates that the magnetic prop- 
erty of the model in high temperature regime is dominated by Curie's law x ^ 1/T. 



C. The strong coupling limit 

For u — > oo, K n (k) goes to zero, from Eqs. (41) we have 

e = k 2 - /i. (53) 

where the external field is set to unit. The A;-sector are completely decoupled with A and ji- 
sectors. At arbitrary temperature, the solutions for the r] n , A n are independent of parameters 
A and ji respectively, which gives rise to Eq. (49) The free energy of the system defined by 
Eq. (44) can be solved by integration by part, 

F/L^O-lll^^] (54, 

where the external field is set to zero. 

We are not able to deduce the specific heat directly from the free energy obtained above 
because the chemical potential is a function of temperature. From Eqs. (34), the density of 
charge rapidity has the form 

P= 2^1 + e (* 2 -M)/r (55) 

Integrating the charge density over k space with the condition Eq. (13), we have an explicit 
expression of the chemical potential, 



7T 2 T 2 n-2 



24/ig 

where /x = n 2 D 2 which denotes fj, at zero temperature. Then the free energy becomes 
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(56) 



F/L = fi D 



1 + 



7T 2 T 2 - 

12// 2 ,- 



3tt 



/'o 



3/2 



7T 2 T 2 - 



(57) 



The free energy for the SU(3) invariant spin chain also has a T 2 dependence [22]. 

Since in thermodynamics S = —dF/dT and = TdS/dT, we find the specific heat at 
low temperature is Fermi-liquid like, 

T 



S — C„ — 



6/r 



(58) 



It is the same as the result of one-component case, since for the strong coupling limit the 
color degree of freedom and the charge degree of freedom are decoupled completely, the 
contribution of color degree of freedom to the free energy vanishes. 



VII. CONCLUSIONS 

In this paper, we have solved one dimensional SU(3) bosons with 5-function interaction 
by means of coordinate Bethe-ansatz method. On the basis of Bethe-ansatz equations we 
first discussed the ground state of the Bose system and found that the ground state is a color 
ferromagnetic state which differs from the SU(3) Fermi system greatly. The configuration of 
quantum numbers for the ground state was given explicitly. The low-lying excitations were 
discussed extensively by both analytical and numerical methods. The energy-momentum 
spectra for three type excitations: holon-antiholon, holon-coloron(<r type) and coloron(cx 
type)-coloron(cr type) were plotted for c = 10 and c = 1. We also discussed the dispersion 
relations of four elementary quasi-particles. 

The thermodynamics of the system were studied by using the strategy of TBA. A revised 
version of Gaudin-Takahashi equations were obtained by minimizing the free energy at 
finite temperature. We found the magnetic property of the system at high temperature 
regime is dominated by Curie's law, and for the case of strong coupling the system possess 
the properties of Fermi-liquid like and its specific heat is a linear function of T at low 
temperature. 
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FIGURES 
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FIG. 1. The density of state in fc-space for the ground state. The distribution changes from a 
histogram to a narrow peak gradually for the coupling from strong to weak. The figure is plotted 
for N = L = 41 and c = 10, 1, 0.1, 0.01. 





4 r 




3.5 




3 - 


ates 


2.5 - 






ity of 


2 


lens 


1.5 - 



A 

- 


A 

- 


r- FF= 


~T~1 , — i 



1 

0.5 



-3-2-10123 
quasi-momentum k 

FIG. 2. The density of state in /c-space for the ground state in the presence of one color rapidity 
by choosing J\ = 0. The distribution changes from a histogram to a narrow peak gradually for the 
coupling from strong to weak. The figure is plotted for iV = L = 100 and c = 10, 1, 0.1, 0.01. 
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FIG. 3. The ground state energy E/L versus the coupling constant lnc for different densities 
D = 1.0,0.75,0.5,0.25. 
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FIG. 4. The holon-antiholon excitation spectrum calculated for N = L = 41 and c = 10 (left), 
c = 1 (right). 
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FIG. 5. The holon-coloron(<r type) excitation spectrum calculated for N = L = 21 and c = 10 
(left), c = 1 (right). 
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FIG. 6. The coloron(<r type)-coloron(<r type) excitation spectrum calculated for N = L = 21 
and c = 10 (left), c = 1 (right). 
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FIG. 7. The dispersion relation of antiholon excitation for different coupling constants where 
the curves from bottom to top correspond to c = 1, 10, 20, 40, 80 respectively. Here N = L = AO. 
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FIG. 8. The dispersion relation of holon for different coupling constants where the curves from 
bottom to top correspond to c = 1, 10,20,40,80 respectively. Here N = L = 40. 



2.5 




1 2 3 4 5 6 



momentum 

FIG. 9. The dispersion relation of a-coloron for different coupling constants where the curves 
from top to bottom correspond to c = 1, 10, 20, 40, 80 respectively. Here N = L = 41. 
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FIG. 10. The dispersion relation of w-coloron for different coupling constants where the curves 
from top to bottom correspond to c = 1, 10, 20, 40, 80 respectively. Here N = L = 40. Zero energy 
corresponds to M = N/2 ground state. 
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